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Abstract 



m ■ 
o ■ 

O ' We discuss the theory of inclusive production of heavy quarkonium, the com- 

parison between theory and experiment, and the expected nuclear effects in 
cold nuclear targets. We also present predictions, based on Nonrelativistic 
QCD (NRQCD), for quarkonium production cross sections in the LHC en- 
ergy range. We find that nuclear effects in the production cross sections are 
largely independent of the sizes of the NRQCD matrix elements, the charmo- 
nium state that is produced and, in the color-octet case, largely independent of 
' the partonic subprocess that produces the heavy quark-antiquark pair. 

> : 
^t- . 

m 1. A Review of NRQCD 

O ■ 

ITj I 11 The NRQCD Factorization Method 

In both heavy-quarkonium decays and hard-scattering production, large energy-momentum scales ap- 
pear. The heavy-quark mass tuq is much larger than Aqcd, and, in the case of production, the trans- 
verse momentum px can be much larger than Aqcd as well. Thus, the associated values of a s are much 
Q_i! less than one: a s (m c ) 0.25 and a s (mi,) 0.18. Therefore, one might hope that it would be possi- 

O ■ ble to calculate the rates for heavy quarkonium production and decay accurately in perturbation theory. 

However, there are clearly low-momentum, nonperturbative effects associated with the dynamics of the 
quarkonium bound state that invalidate the direct application of perturbation theory. 
►> ■ 

• rH , In order to make use of perturbative methods, one must first separate the short-distance/high- 

momentum, perturbative effects from the long-distance/low-momentum, nonperturbative effects — a pro- 
cess which is known as "factorization." One convenient way to carry out this separation is through the 
use of the effective field theory Nonrelativistic QCD (NRQCD) [QEIEI. NRQCD reproduces full QCD 
accurately at momentum scales of order itiqv and smaller, where v is heavy-quark velocity in the bound 
state in the center-of-mass (CM) frame, with v 2 « 0.3 for charmonium and » 2 « 0.1 for bottomonium. 
Virtual processes involving momentum scales of order itiq and larger can affect the lower-momentum 
processes. Their effects are taken into account through the short-distance coefficients of the operators 
that appear in the NRQCD action. 

Because QQ production occurs at momentum scales of order itlq or larger, it manifests itself 
in NRQCD through contact interactions. As a result, the quarkonium production cross section can be 
written as a sum of the products of NRQCD matrix elements and short-distance coefficients: 

^) = E^(0|^(A)|0). (1) 

n rn% 

Here, H is the quarkonium state, A is the ultraviolet cutoff of the effective theory, the F n are short- 
distance coefficients, and the C„ are four-fermion operators, whose mass dimensions are d n . A formula 
similar to Eq. Q exists for the inclusive quarkonium annihilation rate 0. 

* To appear in the CERN Yellow Report for the 2001-2002 CERN workshop Hard Probes in Heavy Ion Collisions at the 
LHC. 



The short-distance coefficients F n (A) are essentially the process-dependent partonic cross sections 
to make a QQ pair. The QQ pair can be produced in a color-singlet state or in a color-octet state. 
The short-distance coefficients are determined by matching the square of the production amplitude in 
NRQCD to full QCD. Because the QQ production scale is of order mg or greater, this matching can be 
carried out in perturbation theory. 

The four-fermion operators in Eq. Q create a QQ pair, project it onto an intermediate state that 
consists of a heavy quarkonium plus anything, and then annihilate the QQ pair. The vacuum matrix 
element of such an operator is the probability for a QQ pair to form a quarkonium plus anything. These 
matrix elements are somewhat analogous to parton fragmentation functions. They contain all of the 
nonperturbative physics that is associated with evolution of the QQ pair into a quarkonium state. 

Both color-singlet and color-octet four-fermion operators appear in Eq. ©. They correspond, 
respectively, to the evolution of a QQ pair in a relative color-singlet state or a relative color-octet state 
into a color-singlet quarkonium. If we drop all of the color-octet contributions in Eq. Q, then we have 
the color-singlet model [4|. In contrast, NRQCD is not a model, but a rigorous consequence of QCD in 
the limit AqcdMq -> 0. 

The NRQCD decay matrix elements can be calculated in lattice simulations 13 16) or determined 
from phenomenology. However, at present, the production matrix elements must be obtained phe- 
nomenologically, as it is not yet known how to formulate the calculation of production matrix elements 
in lattice simulations. In general, the production matrix elements are different from the decay matrix 
elements. However, in the color-singlet case, the production and decay matrix elements can be related 
through the vacuum-saturation approximation, up to corrections of relative order v 4 0. 

An important property of the matrix elements, which greatly increases the predictive power of 
NRQCD, is the fact that they are universal, i.e., process independent. NRQCD w-power-counting rules 
organize the sum over operators in Eq. © as an expansion in powers of v. Through a given order in v, 
only a limited number of operator matrix elements contribute. Furthermore, at leading order in v, there 
are simplifying relations between operator matrix elements, such as the heavy-quark spin symmetry 
and the vacuum-saturation approximation [3], that reduce the number of independent phenomenological 
parameters. In contrast, the CEM ignores the hierarchy of matrix elements in the v expansion. 

The proof of the factorization formula Q relies both on NRQCD and on the all-orders perturbative 
machinery for proving hard-scattering factorization. A detailed proof does not yet exist, but work is 
in progress 0. Collections to the hard-scattering part of the factorization are thought to be of order 
(m,Qv/pT) 2 , not (tuq/pt) 2 , in the unpolarized case and of order m,Qv/pT, not tuq/pt, in the polarized 
case. It is not known if there is a factorization formula at low px or for the pj-integrated cross section. 
The presence of soft gluons in the quarkonium binding process makes the application of the standard 
factorization techniques problematic at low pr. 

In the decay case, the color-octet matrix elements can be interpreted as the probability to find the 
quarkonium in a Fock state consisting of a QQ pair plus some gluons. It is a common misconception that 
color-octet production proceeds, like color-octet decay, through a higher Fock state. However, in color- 
octet production, the gluons that neutralize the color are in the final state, not the initial state. There is 
a higher-Fock-state process, but it requires the production of gluons that are nearly collinear to the QQ 
pair, and it is, therefore, suppressed by additional powers of v. 

In practical theoretical calculations of the quarkonium production and decay rates, a number of 
significant uncertainties arise. In many instances, the series in a s and in v of Eq. © converge slowly, 
and the uncertainties from their truncation are large — sometimes of order 100%. In addition, the matrix 
elements are often poorly determined, either from phenomenology or lattice measurements, and the im- 
portant linear combinations of matrix elements vary from process to process, making tests of universality 
difficult. There are also large uncertainties in the heavy-quark masses (approximately 10% for m c and 
5% for nib, for the mass ranges used in the calculations) that can be very significant for quarkonium rates 
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Fig. 1: J/tp cross section as a function of pr- The data points are from the CDF measurement [8]. The solid curve is the 
NRQCD-factorization fit to the data given in Ref. |9 1. The other curves give various contributions to the NRQCD-factorization 
fit. From Ref. |9). 

proportional to a large power of the mass. 

1.2 Experimental Tests of NRQCD Factorization 

Here, we give a brief review of some of the successes of NRQCD, as well as some of the open questions. 
We concentrate on hadroproduction results for both unpolarized and polarized production. We also 
discuss briefly some recent two-photon, e + e~, and photoproduction results. 

Using the NRQCD-factorization approach, one can obtain a good fit to the high-py CDF data [ 8 1, 
while the color-singlet model under predicts the data by more than an order of magnitude. (See Fig.^) 
The pt dependence of the unpolarized Tevatron charmonium data has been studied under a number of 
model assumptions, including LO collinear factorization, parton-shower radiation, /ct smearing, and 
factorization. (See Ref. |9| for a review.) 

Several uncertainties in the theoretical predictions affect the extraction of the NRQCD charmoni- 
um-production matrix elements from the data. There are large uncertainties in the theoretical predictions 
that arise from the choices of the factorization scale, the renormalization scale, and the parton distribu- 
tions. The extracted values of the octet matrix elements are very sensitive to the small-pr behavior of 
the cross section and this, in turn, leads to a sensitivity to the behavior of the small- x gluon distribution. 
Furthermore, the effects of multiple soft-gluon emission are important, and their omission in the fixed- 
order perturbative calculations leads to overestimates of the matrix elements. Effects of higher-order 
corrections in a s are a further uncertainty in the theoretical predictions. Similar theoretical uncertainties 
arise in the extraction of the NRQCD production matrix elements for the T ffOll states, but, owing to 
large statistical uncertainties, they are less significant for the fits than in the charmonium case. 

At large px (pr ~ 4m c for the J/ tp) the dominant quarkonium-production mechanism is gluon 
fragmentation into a QQ pair in a 3 Si color-octet state. The fragmenting gluon is nearly on mass shell 
and is, therefore, transversely polarized. Furthermore, the velocity-scaling rules predict that the color- 
octet QQ state retains its transverse polarization as it evolves into S-wave quarkonium [ 1 1 ], up to cor- 





Fig. 2: Left-hand side: J/tp polarization at the Tevatron. The band is the total NRQCD-factorization prediction. The other 
curves give the contributions from feeddown from higher charmonium states. Right-hand side: ip' polarization at the Tevatron. 
The bands give various NRQCD-factorization predictions. The data points are from the CDF measurement 1 17 1. From Ref. 1 15 1. 

rections of relative order v 2 . Radiative corrections, color-singlet production, and feeddown from higher 
states can dilute the quarkonium polarization fT"2l IT3l 1741 1731 IT6ll, Despite this dilution, a substantial 
polarization is expected at large px- Its detection would be a "smoking gun" for the presence of color- 
octet production. In contrast, the color-evaporation model predicts no quarkonium polarization. The 
CDF measurement of the J/ip and ip' polarization as a function of px IT71 is shown in Fig.|2j along with 
the NRQCD factorization prediction JT3lfl4l[T5ll . The analysis of ip' polarization is simpler than for the 
J/ip, since feeddown does not play a role. However, the statistics are not as good for the ip'. The degree 
of polarization is a = (1 — 3£) /(l + £), where £ is the fraction of events with longitudinal polarization. 
a = 1 corresponds to 100% transverse polarization, and a = — 1 corresponds to 100% longitudinal 
polarization. The observed polarization is in relatively good agreement with the prediction, except in the 
highest pt bin, although the prediction of increasing polarization with increasing pt is not in evidence. 

Because the polarization depends on a ratio of matrix elements, some of the theoretical uncertain- 
ties are reduced compared with those in the production cross section, and, so, the polarization is probably 
not strongly affected by multiple soft-gluon emission or K factors. Contributions of higher order in a s 
could conceivably change the rates for the various spin states by a factor of two. Therefore, it is impor- 
tant to carry out the NLO calculation, which involves significant computational difficulties. It is known 
that order- corrections to partem fragmentation into quarkonium can be quite large 1181 . If spin-flip 
corrections to the NRQCD matrix elements, which are nominally suppressed by powers of v, are also 
large, perhaps because the velocity-scaling rules need to be modified, then spin-flip contributions could 
significantly dilute the J/ip polarization. Nevertheless, in the context of NRQCD, it is difficult to see 
how there could not be substantial charmonium polarization for pt > 4:m c . 

Compared to the J/ ^-polarization prediction, the T-polarization prediction has smaller w-expan- 
sion uncertainties. However, because of the higher T mass, it is necessary to go to higher px to insure that 
fragmentation dominates and that there is substantial polarization. Unfortunately, the current Tevatron 
data run out of statistics in this high-pr region. CDF finds that a = —0.12 ± 0.22 for 8 < px < 20 GeV 
[19 1, which is consistent with both the NRQCD-factorization prediction [20 1 and the zero-polarization 
prediction of the CEM. There are also discrepancies between the polarizations observed in fixed-target 
experiments and the NRQCD predictions. 

Calculations of inclusive J/ip and T production in 77 collisions |2"T1 l22l have been compared 
with LEP data |23]|24l|25l. Both the J/ip and T measurements favor the NRQCD predictions over those 



of the color-singlet model. 

Belle [26] and BaBar [27| have also measured the J '/ip total cross sections in e + e _ — > J/ipX. 
The results of the two experiments are incompatible with each other, but they both seem to favor NRQCD 
over the color-singlet model. A surprising new result from Belle [28] is that most of the produced J/i/j's 
are accompanied by an additional ccpair: a(e + e~ — ► J/ip cc)/a(e + e~ — > J/xfj X) = 0.591q'J| ±0.12. 
Perturbative QCD plus the color-singlet model predict that this ratio should be about 0.1 K9T There 
seems to be a major discrepancy between theory and experiment. However, the order-a;? calculation 
lacks color-octet contributions, including those that produce J/ tp cc. Although these contributions are 
suppressed by v 4 rj 0.1, it is possible that the short-distance coefficients are large. In other results, 
the angular distributions favor NRQCD, but the polarization measurements show no evidence of the 
transverse polarization that would be expected in color-octet production. However, the center-of-mass 
momentum is rather small, and, hence, one would not expect the polarization to be large. 

Quarkonium production has also been measured in inelastic photoproduction l30l [3 1 1 and deep- 
inelastic scattering (DIS) 1321 133 1 at HERA. The NRQCD calculation deviates from the data near large 
photon-momentum fractions, owing to the large LO color-octet contribution. The NLO color-singlet 
result agrees with the data over all momentum fractions, as well as with the data as a function of px- 
See Ref. [9 1 for a more complete review. In the case of deep-inelastic scattering, the Q 2 and px depen- 
dences are in agreement with NRQCD, but the results are more ambiguous for the dependence on the 
longitudinal momentum fraction. 



1.3 Quarkonium Production in Nuclear Matter 

The existing factorization "theorems" for quarkonium production in hadronic collisions are for cold 
hadronic matter. These theorems predict that nuclear matter is "transparent" for J/ip production at large 
PT- That is, at large pr, all of the nuclear effects are contained in the nuclear parton distributions. The 
corrections to this transparency are of order {rriQv/pT) 2 for unpolarized cross sections and of order 
ijiqv/pt for polarized cross sections. 

The effects of transverse-momentum kicks from multiple elastic collisions between active partons 
and spectators in the nucleons are among those effects that are suppressed by (itiqv/pt) 2 - Nevertheless, 
these multiple-scattering effects can be important because the production cross section falls steeply with 
Pt and because the number of scatterings grows linearly with the path length through nuclear matter. 
Such elastic interactions can be expressed in terms of eikonal interactions [34] or higher-twist matrix 
elements [35 1. 

Inelastic scattering of quarkonium by nuclear matter is also an effect of higher order in (rnQv/px) 2 - 
However, it can become dominant when the amount of nuclear matter that is traversed by the quarkonium 
is sufficiently large. Factorization breaks down when 

> min(z Q ,Zq)P£ 
~ M A k T (tot) ' 

where L is the length of the quarkonium path in the nucleus, Ma is the mass of the nucleus, z is the 
parton longitudinal momentum fraction, Pjj is the momentum of the quarkonium in the parton CM 
frame, and fey (tot) is the accumulated transverse momentum "kick" from passage through the nuclear 
matter. This condition for the breakdown of factorization is similar to "target-length condition" in Drell- 
Yan production 1361 l37l . Such a breakdown is observed in the Cronin effect at low pt and in Drell-Yan 
production at low Q 2 , where the cross section is proportional to A a , and a < 1. 

It is possible that multiple-scattering effects may be larger for color-octet production than for 
color-singlet production. In the case of color-octet production, the pre-quarkonium QQ system carries 
a nonzero color charge and, therefore, has a larger amplitude to exchange soft gluons with spectator 
partons. 



Table 1: Matrix elements for charmonium production. Note that here (£>f } = (Of ( - nS) ( i S 1 )} for J/ip and if)', but (Of) = 
(Oj c0 ( 3 Po)} for XcO- Uncertainties are statistical only. From Ref. 1151 . 
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1.3 ±0.1 GeV 3 


(4.4 ±0.7) x 10~ 3 


GeV 3 
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(8.7 ± 0.9) x 10~ 2 GeV 3 


if 


0.65 ± 0.06 GeV 3 


(4.2 ± 1.0) x 10~ 3 


GeV 3 


3.5 


(1.3 ±0.5) x 10~ 2 GeV 3 


XcO 


(8.9 ± 1.3) x 10~ 2 GeV 5 


(2.3 ±0.3) x 10~ 3 


GeV 3 







At present, there is no complete, rigorous theory to account for all of the effects of multiple 
scattering and we must resort to "QCD-inspired" models. A reasonable requirement for models is that 
they be constructed so that they are compatible with the factorization result in the large-p^ limit. Many 
models treat interactions of the pre-quarkonium with the nucleus as on-shell (Glauber) scattering. This 
assumption should be examined carefully, as on-shell scattering is known, from the factorization proofs, 
not to be a valid approximation in leading order in (itiqv/pt) 2 - 

2. NRQCD Predictions for the LHC 

In this section, we shall use the formalism of NRQCD to give predictions for quarkonium production 
in the LHC energy range. We rewrite the cross section in Eq. Q for the inclusive production of a 
charmonium state H as follows: 

a(H) = Y,° m)n (0%) , (3) 

n 

where <7^Q)n = F n (A)/m d Q n ' 4 , {0%) = (0|Of |0), and n runs over all the color and angular mo- 
mentum states of the QQ pair. The cross sections <t(W)» can be calculated in perturbative QCD. All 
dependence on the final state H is contained in the nonperturbative NRQCD matrix elements (Of). 

The most important matrix elements for J/ip = ip(lS) and ifi' = ip(2S) production can be re- 
duced to the color-singlet parameter (of^ nS \ 3 Si)} and the three color-octet parameters (Of ^ nS \ 3 Si)}, 

(0^ nS \ l So)} , and (0^ nS \ 3 Pq)} ■ Two of the three color-octet matrix elements only appear in the 
linear combination 

Mt {nS) = (k/m 2 c )(ot {nS \ 3 P )} + (Ot^CSo)} . (4) 

The value of k is sensitive to the px dependence of the fit. At the Tevatron, i«3. Fits to fixed-target total 
cross sections give larger values, k « 6-7 [38]. The most important matrix elements for Xcj production 
can be reduced to a color-singlet parameter (O* c0 ( 3 Pq)} and a single color-octet parameter (Og c0 ( 3 Si)). 
These matrix elements are sufficient to calculate the prompt J/ip cross section to leading order in a s and 
to order v A relative to the color-singlet contribution. 

In pp collisions, different partonic processes for QQ production dominate in different pt ranges. 
If px is of order itiq, fusion processes dominate, and, so, the QQ pair is produced in the hard-scattering 
process. These contributions can be written in the form 

a Fu (H) = Y.J dx ± dx * fi/A(x l ^ 2 )f j/B (x 2 ,p J 2 )4f )n (Of) , (5) 




Fig. 3: Differential cross sections per nucleon multiplied by leptonic branching fractions for prompt J/ip (upper curves), J/ip 
from Xc (middle curves), and prompt ip(2S) (lower curves) in pp, pPb, and Pb+Pb collisions at t/s = 5.5 TeV. The EKS98 
parametrization |48 49 1 is employed for pPb and Pb+Pb collisions. 

where A and B are the incoming hadrons or nuclei. In Eq. ©, we include the parton processes ij — > 

QQ X, where ij = gg, qq, qg and qg, and q = u,d, s. The relevant partonic cross sections a^^' n are 
given in Refs. 1391 14011 . 

For pt 3> rriQ, the dominant partonic process is gluon fragmentation through the color-octet 3 Si 
channel. This contribution can be expressed as 

a Fv (H) = E/ da; i dx * dzf i/A { Xl ,^)f j/B {x 2 ^ 2 )afpf^ s ^{z^l T ) (C%) , (6) 

where D? Q ^\z , //| r ) is the fragmentation function for a gluon fragmenting into a QQ pair, P/z is 
the momentum of the fragmenting gluon, P is the momentum of the QQ pair, and /ip r is the fragmen- 
tation scale. The fragmentation process scales as da /dp T ~ IS El- The fragmentation process 
is actually included in the fusion processes of Eq. ©. In the limit px 3> tuq, the fusion processes that 
proceed through g* — > (QQ)s( 3 Si) are well-approximated by the expression ©. At large pr, one can 
evolve the fragmentation function in the scale ^p r , thereby resumming large logarithms of p^/rriQ. Such 
a procedure leads to a smaller short-distance factor fTHl and a more accurate prediction at large p? than 
would be obtained by using the fusion cross section (J5J). However, in our calculations, we employ the 
fusion cross section ©, which leads to systematic over-estimation of the cross section at large px- 

In order to predict the cross section for prompt J/ip production (including Xc and ip' feeddown) 
at the LHC, we need the values of the NRQCD matrix elements. There have been several previous 

extractions of the color-octet matrix elements 031[r5l|39l|40l|43l|Sl|45) from the CDF Xc and ^ 
Pt distributions IHl l4*r51l . We use the matrix elements given in Ref. |T5J, which are shown in Tabled Our 
calculations are based on the MRST LO parton distributions \41\. In calculating the cross section per 
nucleon for prompt J/ip production in pA or AA collisions, we take f^i A = f^pRf. We employ the 
EKS98 parametrization [48 49] for the nuclear shadowing ratio Rf. We evolve a s at one-loop accuracy, 
and we set fi = (4m;? + p T ) l l 2 and m c = 1.5 GeV. 

There are several sources of uncertainty in our predictions for the cross sections. There are large 




Fig. 4: The pr dependence of Rab [Eq. {j]. (a) We compare the results for pPb and Pb+Pb collisions. (The pr dependence 
is stronger in the Pb+Pb result.) (b) We show the dependence of Rab on the various production channels in Pb+Pb collisions 
at yfs = 5.5 TeV. 

uncertainties in the NRQCD matrix elements themselves. The errors shown in Table [2 are statistical 
only. There are additional large uncertainties in the matrix elements that arise from truncations of the 
series in a s and v in the theoretical expressions that are used to extract the matrix elements. The matrix 
elements (Os( 1 S'o)) and (08( 3 -Po)) are fixed by the data only in the linear combination Mj? . In the 
present calculation, we take {0%(}Sq)) = iMf and (Os( 3 Pq)) /ml = (1 — x)M^ /k, use the values 
of k given in Table and choose x = 1/2. Variation of x between and 1 affects the cross sections 
at low px by amounts on the order of 5%. There are additional uncertainties in the predicted cross 
sections that arise from the choices of the parton distributions, the charm-quark mass m c , and the scale 
fi. Because they affect the matrix-element fits, these uncertainties are highly correlated with those of the 
matrix elements. We have not tried to estimate their effects on the predicted cross sections. 

In Fig. |3j we show the px distributions per nucleon multiplied by the dilepton branching fractions 
for prompt J/tp (upper curves), J/ip from Xc decays (middle curves), and prompt ifi' (lower curves) at 
= 5.5 TeV. For pPb and Pb+Pb collisions, we use the EKS98 parametrization I48l l49ll to account for 
the effect of nuclear shadowing. The pp, pPb, and Pb+Pb results essentially lie on top of each other in 
Fig. EJ owing to the many decades covered in the plot. 

In order to display small differences between the distributions, we define the function Rab- 

daAB/dpT - da P p/dpT 

Rab{Pt) = 1 -n • (/) 

da pp /dp T 

In Fig. we present Rab as a function of px- As is shown in Fig.|4la), nuclear shadowing increases 
the cross section at large px and decreases it at small px- The deviation of the Pb+Pb cross section from 
the pp cross section is twice as large as that seen in the case of pPb collisions. In order to investigate the 
dependence of the shadowing effect on the short-distance cross sections that arise in hadroproduction of 
S-wave charmonium states in Pb+Pb collisions, we plot Rab for all channels separately. [See Fig.|4lb).] 
Even though the px dependence the contribution to the cross section of the color-octet 3 Si channel is 
quite different from those of the color-octet 1 So and 3 Pj channels, all three channels show the same 
nuclear effect. The only channel that shows a slightly different behavior is the color-singlet channel, 
which gives a negligible contribution to the cross section. While the differential cross sections in Fig.|3] 



Table 2: Inclusive color-singlet matrix elements for bottomonium production. The errors on the 3 Si matrix elements come 
from estimates of the T(nS) decay rate to lepton pairs. The errors on the 3 Pj states come from an average over potential- 
model estimates. The inclusive matrix elements are a linear combination of branching ratios, as in Eq. |9j. The S-state matrix 
elements are in units of GeV 3 while the Pj-state matrix elements are in units of GeV 5 . From Ref. 1101 . 
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<0f( 3 Si)> mc 


(Of( 3 P ))inc 


i(Of( 3 P!)) inc 


i(Of( 3 P 2 ))inc 


T(3S) 


4.3 ±0.9 











T(25) 


5.0 ±0.7 


0.12 ±0.06 


0.55 ±0.15 


0.42 ±0.10 


T(1S) 


12.8 ±1.6 


< 0.2 


1.23 ±0.25 


0.84 ±0.15 



are strongly dependent on the nonperturbative NRQCD matrix elements, Rab is almost independent of 
the matrix elements, making it a good observable for studying nuclear shadowing at the LHC. 

The T rates are somewhat more difficult to calculate because of the many feeddown contributions. 
The matrix elements are also not particularly well known. Since it is unlikely that all the different 
contributions can be disentangled, we follow the approach of Ref. [ 10 1 and compute the inclusive T(nS) 
production cross section 

da(T(nS)) mc = da^)^ 3 ^)(0 1 T(nS) ( 3 5 1 )) inc + ^^( 65 )^ 3 ^)(0 1 T(nS) ( 3 P J ))^ 

J 

+daW^(Oj inS) eS 1 )) mc + da^^iO^eSo))^ 

+ (E( 2J + l)^ )8(3Pj) ) (^^(^O))!^ , (8) 

where the last term makes use of heavy-quark spin symmetry to relate all of the octet 3 Pj matrix elements 
to the octet 3 Pq matrix element. The "inclusive" matrix elements are defined by 

(Oj inS \n)) iac =J2B H ^ r{nS) (Of(n)) , (9) 
H 

where i = 1 or 8 for singlet or octet, respectively. The sum over H includes the T(nS) as well as 
all higher states that can decay to T(nS). The branching ratio for H — > H' decays is Bh^h' with 
Bh^h = 1- Only XbO-P) an d X&(2-P) decays are included; the possibility of feeddown from the as-yet 
unobserved Xfe(3P) states is neglected. In the linear combination M^ nS \ the color-octet matrix element 

from the 3 P state is neglected, and, so, M^ {nS) = (Oj {nS) i 1 S )) inc . We use m b = 4.77 GeV and the 
MRST LO partem distributions. The values of the inclusive color-singlet matrix elements are given in 
Table |2j and the values of the inclusive color-octet matrix elements, from Ref. lHOl . are given in Table |3] 

In Fig. 13 we show the pt distributions per nucleon multiplied by the dilepton branching fractions 
for the 3 T S states at = 5.5 TeV. The feeddown contributions are included as in Eq. ©. For pPb 
and Pb+Pb collisions, we use the EKS98 parametrization l48ll49ll in order to account for the effects of 
nuclear shadowing. The pp, pPb, and Pb+Pb results lie essentially on top of each other in Fig. |5] 

The unusual relative behavior of the T(2S) and T(35) states at both low and high pr is due to the 
fact that the bottomonium matrix elements are not very well determined. For pt < 10 GeV, the T(2S) 
cross section drops below the T(35) cross section because the T(25) has a large negative color-octet 
matrix element. (See Table |3]) The short-distance coefficients multiplying M^ nS ^ are significant at 
low pt- Thus, there is a large cancellation between the octet 3 Si matrix element and M^ nS \ which 



Table 3: Inclusive color-octet matrix elements for bottomonium production. The matrix elements were fit using the MRSTLO 
parton distributions. The first set of error bars is from \ 2 fits to the T pr distributions in the region pr > 8 GeV. The second 
set is associated with the variation of the scales and corresponds to multiplying fi = yj m% + p\, by 2 (upper error) and 0.5 
(lower error). The matrix elements are in units of 10~ 2 GeV 3 . From Ref. 1 10 1 . 



H 


(Of(^l))inc (O£eS )) inc -|(Of(3p )) inc 


T(35) 
T(2S) 
T(1S) 


3.7±1.7±^ 7.5 ± 4.9±|;| 
19.6 ±6.3+^| -8.7 ±ll.l"f;| 
11.7±3.0il;^ 18.1±7.2^ 1 1 4 




Fig. 5: Differential cross sections per nucleon multiplied by leptonic branching fractions for inclusive T(IS') (upper curves), 
T(2S) (middle curves), and prompt T(35) (lower curves) in pp, pPb, and Pb+Pb collisions at yfs = 5.5 TeV. The EKS98 
parametrization |48 49 1 is employed for pPb and Pb+Pb collisions. 
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Pt (GeV) p T (GeV) 

Fig. 6: The pr dependence of Rab [Eq. 0] for T production, (a) We compare the results for pPb and Pb+Pb collisions. (The 
pr dependence is stronger in the Pb+Pb result.) (b) We show the dependence of Rab on the various production channels in 
Pb+Pb collisions at y/s = 5.5 TeV. 

reduces the T(2S) cross section in this region, causing it to drop below the T(35) cross section at low 
Pt- At the high-p^ end of the spectrum, the large value of the 3 S\ T(25) color-octet matrix element 
(Tabled causes the T(2S) cross section to approach that of the T(15). In this region, the color-octet 
3 5i contribution dominates the other channels. Its large matrix element gives the T(25) an unreasonably 
large cross section relative to that of the T(15). The T(2S) rate at px ~ 15 GeV is more reasonable 
because the large and positive 3 Si contribution, and the large and negative M k contribution nearly 
cancel each other. 

Better determinations of the T matrix elements are required in order to make more accurate pre- 
dictions of the NRQCD T-production rates at the LHC. As is shown in the T -polarization analysis in 
Ref. [20], some theoretical predictions have quite large uncertainties even at Tevatron energies, owing to 
our poor knowledge of the matrix elements. However, as is shown in Fig. |5J the ratio Rab is still a good 
measure of the effect of shadowing on T production. The ratio is independent of the T state and is quite 
similar to the J/ip ratio in Fig. |4] The shadowing effect in Pb+Pb interactions may be somewhat less for 
the T at pt ~ 5 GeV than for the J/tp, but the difference is small. Note also, from Fig.|6tb), that Rab is 
essentially independent of the matrix elements and is, therefore, largely unaffected by their uncertainties. 
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